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SUMMARY 

The dimensional stability of the surface of a large reflector antenna is 
important when high gain/low sidelobe performance is desired. If the surface 
is distorted due to thermal or structural reasons, antenna performance can be 
improved through the use of an array feed. The design of the array feed and 
its relation to the surface distortion are examined in this paper. The 
sensitivity of antenna performance to changing surface parameters for fixed 
feed array geometries is also studied. This allows determination of the 
limits of usefulness for feed array compensation. 


*Work performed under NASA Grant NAG3-419. 



1 • Introduction and Overview 

Maintaining the surface accuracy of a large reflector antenna is 
important when high gain/low sidelobe performance is desired. Reflectors 
surface errors are generally classified into two following types: 

(i) The random surface error in the order of several mils due to 
manufacture imperfection [1-3]. 

(ii) The large-scaled surface error in the order of inches due to 
structural or thermal causes [4] [5]. 

It is the large-scaled surface error that is of interest here. Methods for 
reducing this error include the use of better antenna supporting struc- 
tures, better thermal-insensitive material, and mechanically tunable sur- 
faces. A different approach is to compensate the degraded antenna 
performance due to surface distortion by the placement, design, and excita- 
tion of the feed. The latter approach is becoming increasingly attractive 
due to recent advances in monolithic microwave integrated circuits. Feed 
arrays with many elements can be or will be built compactly and at a reason- 
able cost. 

At any rate, surface error reduction methods and surface error compen- 
sation techniques are complementary; each approach can further improve 
antenna performance independently of the benefits of the other. In this 
paper, surface error compensation techniques will be examined. Analysis 
will begin with the assumption that the surface distortion is given; then, 
methods will be investigated to determine the compensating feed array 
geometry and element excitations. The sensitivity of the antenna perfor- 
mance to changing surface distortion parameters, for fixed feed array 
geometries, will then be examined. 
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For a better understanding of reflector surface error compensation, it 
is convenient to think of a reflector antenna in the receive mode. Figure 
la depicts a plane wave incident on a perfect parabolic reflector. In the 
geometric optics limit, all rays will converge to a single point at the 
focus of the parabola, where they are collected by a single feed element. 
Figure lb shows a plane wave incident on a distorted parabolic reflector. 

In this case, the rays will not converge to a single point, but will be 
spread over the focal plane. To provide reflector distortion compensation, 
a feed array must be designed that will capture the energy spread over the 
focal plane and add it constructively in the feed network. Improvements in 
antenna performance due to compensation in the receive mode will produce 
improvements in antenna performance in the transmit mode by reciprocity. 

Before designing a f6ed array to compensate for surface distortion, we 
must first determine what the important properties of the surface distor- 
tion are and how they affect the geometry of the feed array. This will 
determine the limits of usefulness for feed array compensation. In this 
section, we will give a simple approximate formula, showing that the ampli- 
tude and spacial frequency of the distortion determines the size of the 
feed array necessary for compensation. If the spacial frequency and ampli- 
tude of the surface distortion can be predicted through structural analysis 
or measurement, an estimate of the feed array area can then be obtained. 

Referring to Figure lc, let us consider an off-set parabolic reflector 
with focal point at F. The geometry of the undisturbed reflector is 
completely specified by angles: 

0 = off-set angle of the reflector 
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half-angular extent of the reflector surface 


0 

max 


Along the y-direction, the surface distortion is described by two parameters: 


0 = the amplitude of distortion measured radially from the feed, or 

parallel to the line connecting the vertex V to the focal point F. 
N » # of periods in the distortion over the extent of the reflector 

surface 


Then the minimum array size d in the y-direction needed for compensating 
the surface error is approximately given by 


d 


eirNg 
sin 0 


max 


(i + cos 0 off ) 


( 1 ) 


where e » 2.718 tt = 3.141 d and 0 have the same length unit, e.g. 

wavelength. Derivation and more detailed explanation of (1) is given in 
Section 2. 


As an example, assume that 0 = 20° and 0 rt = 45°. Also assume that 

r max 0 

the distortion on the reflector can be approximated by sinusoidal shape 
along y-direction with a period of 1.0 and an amplitude of 0.30 X. Then 
the y-direction of the feed array will need to be at least 


d 


eir(0.30)( 1 .0) 
sin (20°) 


[1 + cos (45°)] = 12.8 


X 


to fully compensate for the surface distortion. If the distortion Is not 
well approximated by a single sinusoid, it can always decompose into a sum 
of sinusoids. The component in the summation with the largest 0N product 
could then be used to get an estimate of necessary feed array size. Note 
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that if the distortion amplitude , or the number of periods over the extent 
of the reflector surface were doubled, the feed array size in that direc- 
tion would have to be doubled for complete compensation. 

For a complete determination of the feed geometry, the feed element 
spacing must also be determined. It will be shown that element spacing is 
related to the focal length and diameter of the reflector. 

Once the feed array geometry is known, the amplitude and phase excita- 
tions of each array element must be determined. Section 3 will describe two 
methods for the determination of element excitations: the conjugate field 
matching technique and the optimum directivity method. The estimated feed 
geometry and the given reflector distortion profile are combined with the 
physical optics approximation to calculate field quantities in both 
methods. Thus, a fairly accurate estimate of antenna performance can be 
obtained, allowing a fine tuning of the feed array geometry for a given 
distortion profile. 

In practice, the reflector distortion profile is usually not known 
exactly, and may even be a function of time or orientation of the antenna. 

It would therefore be necessary to examine how performance, for a given 
antenna design, is affected by changes in surface distortion. Section 4 
utilizes the results of Section 3 to determine how boresight gain is degraded 
as a function of reflector surface distortion amplitude and spacial fre- 
quency for fixed feed array geometries. The results show how sensitive 
antenna directivity is to changes in the distortion profile. 
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2. Determination of Feed Geometry 

We shall now derive formula (1), which estimates the minimum array 

dimension needed for compensating the effect of surface distortion • 

Consider the PEC (perfect electric conductor) reflector shown in Figure 2. 

An incident plane wave traveling in the +2^ direction induces electric 

currents 3^ on the reflector surface, which in turn produce the scattered 

fields (£ , S ) that we wish to determine. This problem can be separated 
s s 

into two parts. The first part concerns the calculation of (? , if ) on the 

s s 

z - z q plane. Here we will assume the reflector is large and smooth enough 
so that geometric optics can be used to determine the scattered fields. In 
the geometric optics limit, a section of a perfect paraboloid will produce 
a truncated incoming spherical wave that will converge at the focal point 
of the paraboloid. If the reflector is distorted, we seek to relate the 
distortion to effective amplitude and phase errors across the z 3 z q plane. 
The second part of the problem is to determine (? g , if g ) for z < z q once the 
scattered fields are known over the z * z^ plane. (£ g , lf g ) can be calcu- 
lated for z < z q by employing the principles of equivalent sources and 
images. It should be observed that in order to apply geometric optics, it 
is assumed that the distortion is smooth enough so that geometric optics 
can be applied. For distortions that have many local maxima and minima 
that are not separated by a large number of wavelengths , this assumption 
breaks down. 

In order to gain an intuitive understanding of how surface distortions 
affect the fields on the z - z q plane, consider the reflector surface 
distortion, AR(0,<j>), that is superimposed on a perfect paraboloid. The 
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distorted reflector surface is then given by 


R"*( 9» <fr) - R (e) + AR (8,$) =■ 2f (1 ~ C ° S - 9 - ) - + AR (9,<fr) (2) 

sin^0 

where R'(0,<|>) = distance from the focal point F to a point on the 

distorted reflector 

R (8) =» distance from F to a point on the perfect paraboloid 
f * focal length of the paraboloid. 

It is shown in Appendix A that the phase error, Ay(c,8,<fr), resulting from 
the distortion AR(0,<fr) has maxima, minima, and saddle points occurring at 
approximately the same values of 0,$ as the maxima, minima and saddle points 
of AR(9,<fr). Furthermore, these extreme values are related by 

r\ 

AT(c,e i> <)) i ) 5—^ AR(9 i ,<t> i ) [1 + cos 0 ± ] (3) 

where (0^,^) * value of 9 and $ where extreme values of AR(0,<j>) occur 
c * constant. 

Although Ay(c,9,$) and AR(9,$) do not have the same shape, they do have 
approximately the same frequency content, and amplification of their 
extreme values is related by [1 + cos 9^] . Note that the variation of AR 
with 9 is assumed large compared with the variation of cos 9 with 0. These 
properties will prove useful in relating the surface distortion with its 
resultant phase distortion. 

The fields across the z * plane (see Figure 3) can be expressed as 
the geometric optics (GO) field multiplied by an amplitude distortion 
weighting function, where the effective phase distortion is added to the 
phase of the GO field. The amplitude distortion function arises from a 
change in the undistorted amplitude distribution caused by the curvature of 
the surface distortion. The scattered electric field can be written for 
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the x component as 


pjkr 0 

E (r ,0 ,<fc )— — - 
sx o o’^o r n 


x e jkr o +jAY(r o >9 o > *o ) 
^x (r o* 9 o**o* A GOx (r o ,0 o’*o* r Q 


z * z 


z = z 


(4) 


jkr 

Where A GOx ( W 


= GO scattered field for perfect paraboloid 


z = z 


A Dx <W*o> 


z = z 


e jAY(9 0 ,Vr o ) 


amplitude distortion weighting function 


phase distortion 


z = z 


A similar expression exists for the y-coraponent . Ay( 9 q , <j> Q ,r Q ) will be 
approximated using the properties of its relation with AR( 0 Q , d> 0 ) * 


The fields for z < z q (Figure 3) can be calculated from the tangential 
component of the electric field on the aperture plane z ■ z q . In Appendix 
B, it is shown that 

A Px ( V W A GOx (r o >9 o ><<> o ) .1kAv(9r>,4>n,r n ) 


E xs ^l^P 



I 


A - f x 2 - fy 2 


a -jk(f xX1 + fyyi) dfx dfy 


(5a) 


E (x. ,y, ) = 
ys 1 1 2-rr 


^y^o * 9 o * ^o^ jkAy( 0 O , <J> 0 » r o) 


A - f x 2 - f y 2 


I 


2 J k(f x x i + f yyi>df x dfy 


(5b) 


8 



\ 


where f 

X 

= sin 0 

o 

COS 6 

o 

f 

* sin 9 

sin <b 

y 

o 

r o 

-► 

r 

* t 

= X X + 

y y + z i 

o 

o 

J o J 0 

-► 

r i 

A 

■ X + 

y i y 


This result was obtained by Rudge et al. [1] using scalar diffraction 
theory. The derivation in Appendix B is a vector formulation using the 
same assumptions as in [1]. 

The x-coraponent of the electric field in the focal plane is the Fourier 
transform of A^ A GQx e divided by the amplitude weighting function 

. A similar result holds for the y-coraponent. For a perfect 
parabolic reflector with a large focal length to diameter (f/d) ratio, 
point source feed, and y-component of electric field equal to zero, 

A Dx e * 1 (no distortion) 


/l-f x 2 —fy 2 


GOx 


1 for the range of f x and fy, ■— > 2 
constant 


and the focal plane distribution reduces to the Fourier transform of the 
aperture function E. Note that for this same reflector in the transmit 
mode, the far-field pattern due to a point source located at the focal 
point will also be the Fourier transform of the aperture function E. This 
is the well-known similarity between the far-field pattern and focal-plane 
distribution. It is interesting to note that this similarity still holds 
for small amounts of surface distortion. 

Geometric optics and equations (5) can be used to calculate the 
spread of power in the focal plane. Once this distribution is known, the 
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feed array shape can be choosen to capture most of the incident power in 
the focal plane. It should be noted that the focal plane power distribu- 
tion can also be calculated by formulating the transmit problem and 
applying reciprocity. 

To assess the effects of surface distortion on the focal plane distribu- 
tion, it would be useful to have a closed-form expression that gives a reason- 
ably accurate value for the spread of power in the focal plane, given the 
properties of the surface distortion profile. To this end, we will use the 
approach developed by Rudge et al. [6], which assumes a separable aperature 
distribution given by 


^x ^o^o^ A GQx (9 o ,(|> o ) e -jkAy( 0 O , $ Q ) 
Jl _ f 2 V 2 

VI t X Ey 


exp f J [ ot x cos(5^i)+a y cos(fz^!z) ] j 

( 6 ) 


where h is the maximum value of f and f in the aperture Z (h * sin 0 ), 

x y max 

and Z has a square shape in f and f . a and a are the amplitudes of 

x y x y 

the phase distortion. Although (5) seems to be an over-simplification, we 
can think of decomposing Ay into its frequency components, such that 


2 -jkA Y _ » * exp p[« xp co * N ip ! C * - + Oy q cos(^Zi^Z)l 

p=l q=l 


} 


(7) 


The Fourier transform of A^ x ^qqx e will equal the convolution of all 

Fourier transfer™ of coa!!!»£i a V cosAs^ ehe Foutler 

transform of A^, and the Fourier transform of A qq x « Generally, distortions 
caused by temperature gradients and gravitational effects will be dominated 
by a small number of spacial frequency components, and the product of (7) 
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can be truncated after a few terms. Substituting the right-hand side of 
(6) into (5a) yields 


yxj.yj) / h e iaxCOS ^~~^ e ^ kf * x ldf x / h e iayCos( Nw1,E 

— h — h 


^ e-J kf yn df 


( 8 ) 


We need to evaluate the integral e i a x cos ( ^ ^ k ^ xX *df x , (the 

— h 

integral over f has a similar form). Evaluating yields 
,N Y Trf 


r h e -a x cos( 2ElJL ) e -jkf x xi df 
-h n J 


sin khx. 

- 2h J.(B..) - + 2h\ (i) n J (a ) 


o Vk x' khx 


1 


’2, a) v 

n -0 


sin[khx.-nN x] sin[khx,+nN x] 

* X ^ XX 


khx,-nN it 
1 x 


khx,-nN it 
1 x 


(9) J 


If the magnitude of the phase distortion a - 0, a sine function is 

obtained. As a increases, the sine function amplitude decreases as a 
X nN x 

function of J Q (a x ), and other sine functions displaced by ± — ^ in the 

direction emerge with amplitudes governed by J ( <* x ) . The emergence of 
these sine functions corresponds to the spreading of power in the focal 
plane as the amplitude of the distortion is increased. Note that the 
larger N x is, the larger the distance is between sine functions for given 
values of a x * Hence, if the distortion has a high spacial frequency content 
(short correlation period), the power will be spread over a larger portion 
of the focal plane than if it had a lower spacial frequency content for the 
same value of a x * We can quantify this observation by determining how 
large n has to be for a given a x so that | J n (ot x ) f is negligible. For n>B>2 
we have 


11 




. (Ox/2)^ 2 

t 'ST ) 


( 10 ) 


Therefore, |j (a ) | goes rapidly to zero for n > 


ea x 


and we can write 


n N ir 3 k x. h 
max x lmax 


ea 

x K T 

— V 


do 


where we assume 


^nmax^ a x^ 


sin[khx, -n N it] sin[khx, +n N tt] 
lmax max x _ lmax max x 


kh -n N ir 
xmax max x 


khx , +n N ir 
lmax max x 


is the last contributing pair of sine functions in the focal plane as n 

increases. We now rearrange (11) to obtain an expression for the width 

d f in the focal plane over which the power will be spread given a and N . 
IX X X 


eci N X 

. « x x 

x * x lmax " 2 sin 9 


( 12 ) 


max 


Since the spacial frequency of the surface distortion profile is close to 
the spacial frequency of its resultant phase distortion, and phase distortion 
amplitude can be approximately related to surface distortion amplitude by 
equation (3), we can write 
eirNg 


d 2 


x " sin 0 


(1 + cos 9^) 


(13) 


max 

Here f$ x is the peak amplitude of the given frequency component of surface 
distortion (in X), and 9^ can be approximated as the angle, in Figure 1, 

between the z-axis and the line connecting the center of the reflector to 
the focal point. (0^ = 0 for prime focus reflectors; 0^ * 0 for offset 
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reflectors.) Equation (13) gives a simple relation between the surface 
distortion properties and how large the feed array needs to be to capture 
the incident power in the focal plane. As we shall see, this formula gives 
a fairly good estimate of array size when we compare the results with a 
more exact numerical analysis using a separable surface distortion. Note 
that if the estimate of d from Equation (13) is too large, the use of com- 
pensating feed array is not feasible. 

With the feed array shape determined, we must now determine the element 
spacing for the complete characterization of the feed geometry. The feed 
element spacing should be as large as possible provided that no grating 
lobes fall inside the angular extent of the reflector as measured from the 
focus. This choice will minimize the number of feed elements without 
significantly affecting antenna performance, as will be shown in the calcu- 
lated results. 
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3. 


Determination of Feed Element Excitations 


With an understanding of how the geometry of the feed array is esti- 
mated from the surface distortion, it is now possible to determine the 
amplitude and phase excitations of the feed elements for compensating 
degraded antenna patterns. In this analysis, it will be convenient to 
think of the reflector antenna in the transmit mode. Let us consider a 
distorted parabolic reflector antenna with a feed array of N elements. The 
complex excitation coefficients, represented by a column vector I, 

I - [Ij, I 2 , ... I N ] T (14) 

are to be determined. To this end, let us first calculate the far field 
(secondary pattern) of the distorted reflector due to the n^ element, 
namely due to the following array excitation 

1=1, and all other I =0 (15) 

m n 

That far field is denoted by 
-jkr 

E (?) [* V 9 ’*> + 5 V 9 **)] (16 > 

Here ( r , 0 , <J>) is the spherical coordinate of the observation point. (it,£) 
are unitary complex vectors describing the (reference, cross) polarizations 
of the antenna. An example of E^ with m = 2 is sketched in Figure 4. The 
column vector 

EO,*) = [E lt E 2 , ... e n ] T (17) 

represents the reference-polarization radiation in direction (0,<J>) of the 
distorted reflector due to each individual element in the feed array. In 
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the present paper, we calculate S(0,<j>) by the standard physical optics 
integration method. 

There are two ways for determining excitation I to improve a degraded 

pattern. The first one is the conjugate field matching method ( CFM) , in 
* 

which we set 


1 ■ ^ (9 ’ 9 beam’ ♦ ' ♦bea»> 


(18) 


Here (0, , ) is the desired main beam direction of the antenna. 

This method has been known for more than ten years, and is widely used in 
practice [5], [10] - [18]. It has been pointed out [15] that, while CFM 
gives good results, it is not an optimum solution in any sense. A truly 
optimum solution which maximizes the antenna directivity in direction 


(9 beam’ <t, beam ) is given by l8] 

I - (IV 1 i*(e - « beall , ♦ - ♦ beam > 


(19) 


Here A is a N x N square matrix describing the coupling between primary 
patterns of the feed elements. The presence of A is due to the fact that 
the total power radiated by the feed array is not equal to the sum of the 
powers radiated by the individual feed elements [19], a statement that 
superposition of power does not hold in general. Now, for the special case 
that spacing between feed elements is large in wavelength, couplings between 

3 = 

element patterns are small. Then A is nearly diagonal and the optimum 
solution in (19) reduces to CFM in (18). 


* Within a proportional constant. 
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The above two methods for determining feed array excitations are orig- 
inally used for compensating degraded secondary patterns arising from off- 
focus feeds. In the next section, we will show that they work just as well 
for compensating the other cause of degraded pattern: reflection surface 

distortion. 
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4. 


Calculated Results 


To outline a method for designing a distortion compensating feed array, 

we will take a specific example of the reflector shown in Figure 5 with 

the distortion shown in Figure 6a and Table 1. Note that parameters (N^, 

Ny) describe ripple numbers along (x,y) directions, respectively. 

Parameters (8,3) describe the corresponding ripple amplitude. First, 
x y 

the focal plane distribution due to a plane wave incident from the bore- 
sight direction is calculated. This can be done using the method of 

Section 2, or applying reciprocity to the solution of the transmitting 

► 

problem. From this analysis, the array of Figure 7a was found to capture 
most of the energy spread over the focal plane. The element spacing was 
choosen such that no grating lobes of the feed array fell on the reflector. 
The optimum directivity method was then applied to determine the element 
weights, and the compensated far-field patterns shown in Figure 8 were 
obtained. These compensated patterns can be compared with the patterns for 
the undistorted reflector using a single feed (Figure 9) and the distorted 
reflector using a single feed (Figure 10). All but 1.7 dB of the 
undistorted reflector gain were recovered with the array feed. 

As noted before, the actual distortion profile is usually not known 
exactly. It would therefore be interesting to see how the antenna gain 
degrades for a fixed feed array geometry when properties of the surface 
distortion are varied. Using the feed geometry in Figure 7 and surface 
distortion of Figure 6a, the degradation of boresight optimum directivity 
with increasing distortion amplitude is shown for three different feed 
array areas in Figure 11a. From these results, it is evident that the feed 
array area increases approximately as the square of distortion 
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amplitude to maintain a given value of gain (i.e., 55 dB) . This is in 
agreement with the analysis of Section 2. Figure lib shows similar results 
for the case when conjugate field matching is used instead of the optimum 
directivity method. Note that the values of directivity are slightly lower 
(the single element case is the same for both the CFM and the optimum methods). 
Figure 12 depicts peak directivity degradation as distortion amplitude 
increases for two arrays of the same area but with a different number of 
elements (i.e., different spacings). Not much improvement is gained by 
using smaller and more closely spaced elements in a given feed array area. 

The three-element array corresponds to the largest element spacing allowed 
without having grating lobes fall on the reflector. Figure 13 shows the 
decrease in peak gain as a function of distortion amplitude for two dif- 
ferent distortion profiles (Figures 6a and 6b). The feed array is the same 
(Figure 7b) in each case. The distortion with higher spacial frequency 
cannot be compensated for as well as the distortion of lower spacial frequency. 
Since the distortions are separable. Equation (1) can be used to estimate the 
maximum value of distortion amplitude at which full compensation can be 
realized. Rearranging Equation (1), we find that this maximum distortion 
amplitude is given by 


d sin (0 ) 

max 

eirN[l + cos 


For the reflector geometry used in the numerical results 

9 = 20 ° 
max 

0 q = 43° 

d = 8 X 


therefore 


18 



0 JL 81 _X 
0 N 

fo.16 X, for N * 1.167 (Distortion 2) 

1 0.24 X, for N = 0.770 (Distortion 1) 

These values are marked by points A and B in Figure 13. Numerically calcu- 
lated gains begin to fall off shortly after these values. This verifies 
the validity of Equation (1). 

In all the cases presented thus far, the feed array amplitude and 
phase weighting are continuous. As a matter of practical importance, it 
would be interesting to know how phase and amplitude quantization affects 
the performance of the antenna. To determine this effect, a plot was made 
(Figure 14) of the degradation in peak boresight directivity as a function 
of quantization. The horizontal axis shows the number of bits of linear 

3 

quantized phase and amplitude. (A 3-bit quantization means 2 steps in 

3 

phase, and 2 steps in amplitude.) This case is for the 25-element feed 

array of Figure 7b and distortion function if 2 with 0 = 8 ■ 0.175 X. Note 

x y 

that a quantization of as few as 3 bits produces a loss of only 0.60 dB 
from the continuous amplitude and phase optimum peak boresight directivity. 
A 5-bit quantization is sufficient for practical purposes. 

From these results, it is apparent that for a given array, the loss of 
gain can be compensated for provided the distortion amplitude and spacial 
frequency remain within certain limits. Therefore, by designing for the 
maximum distortion amplitude and maximum distortion spacial frequency, the 
compensating feed array should work well for the continuum of surface 
distortion profiles up to and including those extremes. 
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5 . 


Conclusion 


The facets of reflector surface distortion compensation were explored 
by first assuming the reflector distortion was given and then designing the 
compensating feed array. We also examined the sensitivity of boresight 
directivity to changing surface distortion parameters for fixed feed-array 
geometries. It was found that feed array compensation is feasible only for 
distortions with low spacial frequency content, such as those induced by 
thermal and gravitational effects. The optimum directivity method for 
determining element excitations pas found to yield slightly better values 
of directivity than those for the conjugate field matching technique. 
However, the CFM technique lends itself to simple realization in hardware, 
and may be very useful in an adaptive implementation. In view of these 
results, distortion compensation using an array feed is indeed a reasonable 
approach to improving antenna performance for large, space-based reflector 
antennas that are not easily accessible for tuning and have time-dependent 
surface distortions. This is particularly true when considering the recent 
advances in MMIC feed array technology. 
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Appendix A: Relation Between Surface Distortion and Resultant 

Phase Distortion 


Let 


R - 2f 


(1 - cosfl) 
sin2e 


r- 2f ° ~, cose) 

sin 2 0 


* distance from F to a point (R,0,$) on the 
perfect parabolic reflector 
+ AR (0,<fr) * distance from F to a point 

(R*,0,+) on the distorted reflector 


The normal to the perfect reflector at the point (R,0,+) has direction 


given by 


* 

V 


. . : 2f (1 - COS0) 2 

r + 9 I 


(A.l ) 


The normal to the distorted reflector at the point (R > ,0,$) has direction 
given by 


> 

V' 


? + 0 


2f (1 - cos0) 2 . 1 3 AR( 0 , ♦ ) , A 1 r 3AR(0,») , 

F gln 3 0 R x 30 1 * R'sinO 1 3* 


(A. 2) 

For values of (0,<fr) where AR(0,$) is a local maximum, minimum or saddle point, 
we have 

3 AR( 0 > — 3 &R( 0 , ♦) m q (A. 3) 

30 3$ 


Therefore, 

> 

V' - - 



(1 - COS0) 2 
sin^0 


(A. 4) 


- r + 0 


It n + 

R 11 R 


(|V - <|V + ...1 


(1 - CO30) 
sin^O 


If « i + * f then V" 2 V, and the normals for the perfect and distorted 
R 

reflectors are approximately equal for values of (0*$) where AR(0»$) 
extreme values. The approximation becomes an exact equality as 0 -► 0. 
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Referring to Figure 15, the difference in path lengths between f PQ f + 
IqU I and | f + | Q"U | is given by 

|pq| + |qu| - |p'Q'| - |q'u| s ar( e, <*>) [1 + cose] (a.5) 

for the values (0,$) corresponding to the extremes of AR(0,$). Since the 
phases at P and P'are equal, then the difference between the phase at u for 
the perfect reflector and the phase at u for the distorted reflector is 
given by 

2 * 

Ay(c,0,$) = ~ AR(0,$) [1 + cos0] (A. 6) 


Now if the sphere r * c is chosen large enough so that none of the reflected 
rays cross before intersecting the sphere, then in the geometric optics 
limit Ay(c,0,$) is at an extreme value on the sphere r * c. This is in 
fact true because a ray is normal to its equiphase surface, and 
7 Ay(r,0,$) « r at r * c. Therefore, 


i 3 Ay(c,9,<») , 0 

c 30 

l 3Ay(c,0,4>) _ 0 

c sin 0 3$ 


(A. 7 ) 
(A. 8) 


Since the phase distribution on r = c due to reflected rays from the per- 
fect parabola is constant, and we are interested only in relative phase on 
the sphere r - c, Ay(c,0,$) is the phase distortion due to the surface 
distortion AR(9,$). 


AR 

+ A typical value for for a sloppy reflector surface is around 

0.001 (assuming AR = 60 mills and R * 60 inches). For these values, 
the change in 8 component will be less than 0.0006 if the reflector 
surface does not extend past 0 * 60°. 

* If AR changes 60 mills for a change in 0 of 25°, the change in 0 
due to the slope of the distortion is around 0.002. 
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Ay(c,0,$) and AR(0t$) have maxima, minima and saddle points occurring 
at approximately the same values of (©,$)• It is therefore reasonable to 
assume that Ay(c,9>$) and AR(9,+) have approximately the same frequency 
content in 9 and Notice also that the maxima, minima and saddle points 
of Ay occur at values independent of r, for r > c. 
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Appendix B: Calculation of Focal Plane Fields 

Consider Figure 2; the induced current on the reflector surface 
produces scattered fields (E g , H g ). can be replaced by the equivalent 
sources J 


n x H and M - E xnonz*z plane to obtain the same 
s s s o 


scattered fields E , H for z < z as produced by the induced current J . . 

s s o i 

♦ * 

By applying the image theory, the fields (E^, H g ) f° r z ^ Z Q can calculated 

-► + 

from the magnetic source M * 2E x n alone* It can be shown [4] (see 

s s 


Figure 3) that 


E(r L ) 


- 7 x 



.' Jkr 10 


2wr 


10 


E(r ) x z da 
o o 


(B.l) 


r 

where E(r^) * electric field on the aperture l 
electric field for z < z 

o 

ds * element of surface area on the z * z plane, 
o o 


o 

ECrj) 


If we approximate E(r^) by a vector weighting function times an 


incoming spherical wave, we have 



jkr 
e J o 


x z ds 


(B .2) 


E (9 ,<k ) x z ds 
s o o o 


sy o o T o 


sy o o T o 


r o ■ 'V 
r l ' I'l 1 


Note from Figure 3 we have, by the law of cosines. 


10 


r + r 
o 


* " 2r i r 0 3in 9 i sin 9 o cos( V*l ) ~ 2r l r o COS 9 1 


cos 0 
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Therefore , 


r, A 2 r,2 

-ifiL . i + _L 

r 2 r.2 


— iy - 2 — [sin sin 0 Q cosC^-^j) + cos 0^ cos 0 q ] 


*• i 

Assuming — « 1 

o 


r , rt 2 - r 2 2 2r,r tsin 0 sin 0 cos($ “♦,) + cos 9. cos 8 1 
10 o 1 o 1 o o T 1 1 o 

r,_ - r 2 r, [sin 0, sin 0 cos(A “♦,) “ cos 0. cos 0 ] 

10 o 1 1 o oi i o 


(B.3) 


x. sin 0 cos 4 + y, sin 0 sin A + z. cos 9 

1 rt ' /\ J I A Q | 


x f + y. f + z f 
lx 1 y 1 z 


where f * sin 0 cos 4 
x o T o 

f * sin 0 sin 4 
y o o 

f - cos 0 
z o 


(B.4) 

(B.5) 

(B.6) 


Substituting (B.3) into (B.2) yields 

-jkCx.f +y,f +z,f ) 


E(r 


, f f -jk(x,f +y. f +z,f ) 

l £ 77 1 y — E , <VW «*] <*- 7 > 


Crossing E (9 ,r ) with z and performing the curl operation yield 
s o o o 




• * /> 


jk(sin 0 o cos 8 o E xa + sin 0 q sin* Q E yg ) z 


-jk(x, f +y, f +z f ) 

+ jk cos 0 E ? + jk cos 0 E 2] — — — ds 

oys oxs rz c 

J o 


(B.8) 


Note that 


ds *o do d4 where * r sin 0 

o o o T o o o o 

z =* r cos 0 
o o o 
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Therefore, 


sin 0 , sin 0 

da m z = d0 da ” r t — d0 di 

o o 3. o o o cos 0„ o T o 

cos 0 o 


(B.9) 


SubstiCuting (B.9) into (B.8) and realizing z^ ■ 0 in the focal plane, we 
have 


« r l> --U 


JJ- 


(tan 0 cos * E + tan 0 sin ^ E ) z 
o^oxs o o ya 


+ E ? + E k] e" Jk(x l f x +y l f y ) sin 0 d0 d* 
ys 7 xa 7 o o T o 


(B.10) 


-ik 

2ir 


JJ - 


111 

/ 1- f - f 

r x y 


E + 
xs 


/ 2 2 

/ 1 - f - f 

Y x y 


E ) z 

ys 


+ E y + E x ] 
ys 7 xs 


1 


e" Jk(x l f x +y lV df df (B.ll) 


111 

/ 1- f - f 

r x y 


x y 
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Distortion 

Function 


Equation 


#1 A z - p x (-0.0000001 574x 4 ) 

+ /J y (12.38 -0.41 38y 
+ 0.0041 00y 2 
- 0.00001 226y 3 ) 

o<p x ,p y <-\.o\ 


#2 AZ = 0 x (1.0- 0.001 953x 2 

+ 0.000000474X 4 ) 
+p y (12.38 -0.4138y 
+ 0.0041 00y 2 
- 0.00001 226y 3 ) 


o<p x ,p y <-\.o\ 


wavelengths @ 30 GHz ) 


Sinusoidal Approximation 


N y tt 

A z * p y cos[ (y- 145.8)] 

ymax 

( x variation neglected ) 

N y = 0.7697 
Vmax - 52 - 8 in - 


U X K 

A Z - P x cos( x) 

x max 

NyJt 

+/3 y cos[ (y-1 45.8)] 

^max 

N x = 1.167 
N y = 0.7697 
x max “ ^max = 52.8 in. 















Ideal Parabolic Surface 



FIGURE 1c. Reflector - Feed Geometry for Equation (1) 




(2=0) 

(b) Equivalent magnetic current M s produces 
scattered fields E s , H s for z less than z 0 . 


FIGURE 2. Geometry for the Problem of Scattering 
from a Reflecting Surface. 
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FIGURE 3. Geometry for the Focal Plane Electric 
Field Calculation. 



Reference 

Polarization 

Radiation 



FIGURE 4. The Reference-Polarization Secondary Pattern in 
direction (e beam ,<t> beam ) is E 2 when only Feed 
Element 2 is Excited with l 2 = 1.0. 



SIDE VIEW END VIEW 



FIGURES. Reflector Geometry Used to Obtain the 

Numerical Results. The feed array is located 

at the z f = 0 plane. The frequency is 30 Ghz (K = 0.3937"). 




y 



(a) distortion function #1 plotted for /? = p « i x 

a y 



(b) distortion function #2 plotted for p x » p y > U 


FIGURE 6. Distortion Functions Superimposed on 
the Perfect Parabolic Reflector Sketched 
in Figure 5. 

(Height dimension in wavelengths @ 30 GHz, base dimension in inches) 
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(a) 



N = 37 Q = 20 

D = 1 .6 X A = 95 X 2 


(b) 


I 


mmnm 

■■via 

nneia 

■ RUB 


N =25 Q = 20 

D = 1 .6 X A = 64 X 2 


I 

* 

I 



N * 9 Q = 20 

D = 1 .6 X A = 23X 2 


(d) 


N = 16 Q= 11.25 

D = 1 .2 X A = 23 X 2 


(e) 


I 

I 

I 

---0--- 

I 

I 


N = 1 Q = 20 

D*1.6X A = 2.5 


Parameter: N = # of elements 

D = element spacing 

Q = element patern Q (E-plane, and H-plane patterns are given by cos Q 0 ) 
A = approximate feed array area 

FIGURE 7. Feed Arrays for the Reflector in Figure 5. 

Each Square Repersents a Feed Element 
Located at the z f = 0 Plane. 


38 






DEGREES 


FIGURE 8. Pattern Cuts for the Distorted Reflector 

with the Array Feed Sketched in Figure 7a. 


RREOEDINO BLANK NOT FILMED 





GRIN 
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DEGREES 


FIGURE 9. Pattern Cuts for the Undistorted Reflector 
with a Single Feed Element Sketched in 
Figure 7e. 
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BORESIGHT DIRECTIVITY (dBi) 



DISTORTION FUNCTION AMPLITUDE p = p - 0 in X 

x y 


FIGURE 11a. Optimum Boresight Directivity as a Function of 
Distortion Amplitude for Three Different 
Array Areas. Distortion function #1 was used 
when calculating each curve. 
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BORESIGHT DIRECTIVITY (dBi) 


I 

i 



DISTORTION FUNCTION AMPLITUDE 0-0 - 0 in X 

x y 


FIGURE 11b. Conjugate Field Matching Boresight Directivity 
as a Function of Distortion Amplitude for 
Three Different Array Areas. Distortion 
function #1 was used when calculating each curve. 



BORESIGHT DIRECTIVITY (dBi) 



DISTORTION FUNCTION AMPLITUDE p =/3 x =/3 y inX 


FIGURE 12. Optimum Boresight Directivity as a Function of 
Distortion Amplitude for Two Different Element 
Spacings ( Equal Array Areas). Distortion function 
#1 was used when calculating both curves. 
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BORESIGHT DIRECTIVITY (dBi) 



DISTORTION FUNCTION AMPLITUDE p = p x = p y \n\ 


FIGURE 13. Optimum Boresight Directivity as a Function of 
Distortion Amplitude for Two Different Distortion 
Functions. The 25 element array of figure 6b was used 
when calculating both curves. 
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Bits Quantization 


FIGURE 14. Effect of weight quantization on peak boresight 
directivity. The weights are quantized linearly 

in phase and amplitude. These results are for the 25 
element feed array (Figure 7b) and distortion function #2 
( p x = p y = 0.1 75 X ) superimposed on the reflector. 
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R'= Distorted z 



FIGURE IS. Relation Between Surface Distortion and 
the Resultant Phase Distortion. 
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